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Abstract 

We investigate the flux flow state in a two-gap superconductor in which two s-wave gaps with 
different amplitudes exist on two separate Fermi surfaces. The flux flow resistivity is obtained 
on the basis of the Bardeen-Stephen relation and the result agrees well with the anomalous field 
dependence of the flow resistivity recently observed in the two-gap superconductor MgB2. Some 
typical properties of the vortex in this system are also discussed. 
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The vortex lines in type II superconductors are subject to the Lorentz force under an external 
current and begin to flow perpendicular to the current and magnetic field when the Lorentz force 
exeeds the pinning force. This is the flux flow state and a finite resistivity arises [1]. The basis of the 
flux flow resistivity is the presence of bound states inside the vortex core. The energy gaps of these 
states are so small that the conductivity in the core is practically normal. Such states were found by 
Caroli et al. by microscopic methods [2]. In both dirty and moderately clean s-wave superconductors, 
the flux flow resistivity pf is proportional to the magnetic field H[l], namely, 
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ti c 2 

where H C 2 is the upper critical field and p n the normal-state resistivity. This is the so-called Bardeen- 
Stephen relation [3]. In the low-temperature region, this relation holds well for almost the entire field 
range of the vortex state. [4] 

The superconductivity of MgB2 has been investigated with keen interest since it has the highest 
transition temperature (T c ~ 39K) among metalic compounds at present and a great number of 
investigations have been carried out [5]. One of the most characteristic features of this superconductor 
is that two s-wave gaps with different amplitudes exist on the two separate Fermi surfaces having 
roughly equal densities of states (DOSs). The two-gap model was proposed on the basis of first- 
principle calculations, [6, 7] and experimental results obtained by point-contact spectroscopy, [8] specific 
heat measurements, [9] and angle-resolved photoemission spectroscopy [10] support the model. 

Recently, the measurement of the flux flow resistivity of MgB2 was reported by Shibata et al.fll] 
Large deviation from the .ff-linear dependence for the flux flow resistivity has been observed, in spite 
of the s-wave pairing symmetry of MgB2[12, 13, 14]. Such anomalous behavior is expected to be 
related to the two-gap structure, but a clear explanation has not yet been proposed. In this study, we 
investigate the flux flow state in the two-gap system and propose a new scenario for the anomalous 
flux flow resistivity. 

First, we investigate the vortex in the two-gap system, before discussing the flux flow state. Let 
\&L and stand for the order parameter for the large energy gap and that for the small energy gap, 
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respectively. We use a Ginzburg-Landau(GL) free energy for the two-gap system in a weak coupling 
approach with a Josephson-type interaction [15], 
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where D = V — ii| A and index _L= We introduce a magnetic field along the z-axis. Let 6>j(r) 
denote the phase of ^(r) (i = L,S). There are two types of vortex, one of which consists of the 
winding of and the other consists of that of 6*s(r). The temperature dependences of the two 

gaps[8, 10] indicate that the mixing effect between the large gap and the small gap should play an 
important role in MgB2- This means that the term proportional to 7 in the GL free energy, eq. (2), 
is not negligible. The term is rewitten as 

- 7 {*J,(r)* s (r) + *L(r)^(r)} = 
-2 7 |* L (r)||* s (r)| cos (0 L (r) - s (r)) , (3) 



and locks the relative phase as 



U £<S: (4) 



Therefore, we may state that (i) the two types of vortex have the same H&, and that (ii) it is favorable 
for their cores to be overlapped energetically, i.e., there is an attractive interaction between the two 
types of vortex. This phase-locking effect could be modified by some boundary effects[16] or in the 
thin film system and this possibility will be discussed elsewhere [17]. 
The coherence length of ^(r) is[l] 



and for |r| >> £l,£S) the amplitudes of order parameters become constant and one obtains a London 
equation 

VxH+(^ + ^)(A-|v„)=0, (6) 

where <I>o = hc/2e, 

mM _n M _ J #s(r) (7 > 0), m 
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A«(T)=(^^|^(T) 1 2 )" /2 , (8) 

and Vl/^(T) is the stationary value of ^i(r) in the homogeneous system at a temperature T. One can 
see from eqs. (6) and (8) that the London penetration depth in this system is 

- / 1 1 \-V2 
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Let us now discuss flux flow resistivity. In the two-gap system, the applied current J is divided 
between the band with the large gap (L-band) and that with the small gap (S-band); 



J = J L + J S . (10) 

The divided current Jj provides the Lorentz driving force to the vortices in the two bands. We assume 
that the two types of vortex have the same velocity because of the presence of the attractive interaction, 
as we mentioned before. The ratio of the two distributed currents is determined to minimize the energy 
dissipation of flux flow, more precisely, the power loss density of the flux flow 

W = W h + W s , (11) 

where 

Wi = irR 2 pfj? (12) 

is the power loss per unit cell of the vortex lattice (with a lattice constant R) per unit length along 
the z-axis[l] and p^ is the flux flow resistivity in the i-band. The field dependence of p^ will be 
discussed later. 

To minimize eq. (11) under the constraint of eq. (10), one obtains 

o {S) 

J, 



o (L) + o {S) 

pf + pf 



Then, 



p (L) 

Js = -(lTT7s) J - ( 13 ) 
Pf + Pf 



and this equation indicates that the total flux flow resistivity in this system is 

1 



" r = w (15) 

which is the parallel connection of the resistivity in the L-band and that in the S-band. 

Let us discuss the field dependence of p^ . Since the two gaps have the s-wave pairing symmetry, 
these two resistivities are considered to obey the Bardeen-Stephen relation, i.e., 

ff = 4>a!? (16) 
for H c i < H < , where p$ is the normal-state resistivity in the i-band, and 

H c2 - (17) 

For MgE>2, the results of the point-contact spectroscopy experiment suggest that the small gap is 
suppressed quicker than the large gap. [8] This implies that £l < £s and < . The value 

coincides with the upper critical field in this system. The kink point that has been observed in 

(s) 

the field dependence of the specific heat [9] corresponds approximately to H c2 , since the presence of 
the kink indicates that the small-gap pairing is almost suppressed [18]. Then, the S-band would quite 
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close to the normal metalic state above H c2 ■ This picture is also supported by the calculation for 
quasiparticle DOS in the vortex state. [19, 20] Therefore, it is possible to approximate 

pf{H)~p& (18) 



for H<§ <H< H^. 



Then, we obtain the field dependence of p™ in eq. (15). By using eqs. (16) and (18), 
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We can see that p 1 ™ is linear with H in the low-field regime, curves convexly in the high field 
regime and is continuously connected to the normal-state resistivity in the two-band system (l/p^ + 
1/^V 1 - 

Let us compare eq. (19) to the experimental results of the flux flow resistivity in MgB 2 [ll]. The 
measurements were performed with the magnetic field parallel to the c-axis of the crystal and also 
with the field in the a6-plane[ll]. The kink point that has been observed in the field dependence of 
the specific heat [9] suggests that H^/H^ = 0.02 for the H \\ ab case and H^f/H^ = 0.1 for the 
H || c case [21]. Unfortunately, there is no definite information for the ratio pffl / Pn^ at present, and 
we use it as a fitting parameter and choose pffl / Pn^ = 0.26 for the H \\ ab case and pffl / Pn^ = 0.5 
for the H \\ c case. As is shown in Fig. 1, it should be emphasized that both cases of the experimental 
results [11] are consistent with the convex- type behavior of eq. (19 ), particularly for the H \\ c case. 
The robustness of the coincidence in the H \\ c case is supported by the fact that the 20% change in 
the parameter pffl / Pn^ around the value 0.5 causes at most a 10% change in p l ™° . The difference 
between our calculated results and the experimental results in the high- field region of the H || ab case 
is considered to be related to complications of the analysis of experimental data, since the vortices 
driven by the microwave field [11] move in the ac-plane and the anisotropy of the crystal causes the 
current-direction dependence of the flux flow resistivity. The experimental data for the H || ab case 
correspond to the averaged value of this anisotropy. 

In summary, we have considered a two-gap superconductor such as MgB2-[5] In this system, there 
are two bands that cross the Fermi level in the metalic phase and two s-wave gaps with different 
amplitudes that arise on these two Fermi surfaces. There are two types of vortex, one of which 
consists of the phase winding of the order parameter for the large gap and the other consists of that 
for the small gap. These two types of vortex have the same lower critical field and attract each 
other. We have examined the flux flow resistivity in this system. The two types of vortex have the 
same velocity because of the presence of the attractive interaction. The flux flow resistivity in the 
two bands is considered to obey the Bardeen-Stephen relation, since the two gaps have the s-wave 
pairing symmetry. We have obtained that, to minimize the power loss (energy dissipation per unit 
time) caused by the flux flow, the total flux flow resistivity is expressed as a parallel connection of 
the resistivity in the band with the large gap and that in the band with the small gap. It should be 
noted that the amplitude of the small gap is suppressed faster than that of the large gap by the 
magnetic field in MgB2[8], and the resistivity in the band with the small gap is approximately equal 
to its nomal state value in the field region where the small gap is almost entirely suppressed. Taking 
this into account, our formula agrees quite well with the anomalous field dependence of the flux flow 
resistivity as observed experimentally[ll]. 

The authors are grateful to J. Akimitsu, H. Fukuyama, N. Furukawa, K. Kubo, K. Ishikawa, 
K. Izawa, Y. Matsuda, S. Soma, P. Thalmeier, and K. Maki for useful discussions and comments. 
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Figure 1: We compare the experimental data of flux flow resistivity in MgE>2 given in ref. 11 and our 
results calculated with eq. (19). Black dots and triangles denote experimental results for the H || ab 
case and H || c case, respectively. The green line denotes eq. (19) with parameters / 'H^ = 0.02 

and ph ] /ph ] = 0.26, and the red line denotes eq. (19) with H^f/H^ = 0.1 and p { n ] / pk ] = 0.5. The 
red line agrees very well with the data for H \\ c. 
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